Excitations of a Bose-Einstein Condensate at Non-Zero Temperature: A Study of 

Zeroth, First, and Second Sound 
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Collective excitations of a dilute Bose gas were probed above and below the Bose-Einstein condensa- 
tion temperature. The temperature dependencies of the frequency and damping rates of condensate 
oscillations indicate significant interactions between the condensate and the thermal cloud. Hydro- 
dynamic oscillations of the thermal cloud were observed, constituting first sound. An antisymmetric 
dipolar oscillation of the thermal cloud and the condensate was studied, representing the bulk flow 
of a superfluid through the normal fluid. The excitations were observed in situ using non-destructive 
imaging techniques. 

PACS numbers: 03.75.Fi, 05.30.Jp, 51.40.+P, 67.90,+z 



Bosc-Einstcin condensation in dilute atomic gases 
has provided a testing ground for well-developed 
many-body theories of quantum fluids. Many studies 
have focused on the collective excitations of such gases, 
which fully describe their dynamics and transport prop- 
erties such as superfluidity. Experiments have studied 
low-lying collective excitations over a range of tempera- 
tures and higher- lying modes H . Zero-temperature 
findings have agreed well with theoretical predictions 
based on a mean-field description of the weak inter- 
atomic interactions |^-|io[]. However, the behavior at non- 
zero temperature, which involves interactions between 
the condensate and the thermal cloud, is not fully un- 
derstood, pointing to the need for new theoretical devel- 
opments. 

The physical nature of collective excitations depends 
on the hierarchy of three length scales: the wavelength 
of the excitation A, the healing length £ which is de- 
termined by the condensate density, and the mean-free 
path / m f p for collisions between the collective excitation 
and other excitations which comprise the thermal cloud. 
The collisionless regime, defined by A -C l m { p , occurs 
at zero temperature and for low densities of the ther- 
mal cloud. There, free-particle excitations are obtained 
at short wavelengths (A <C £) while phonon-like excita- 
tions known as zeroth sound are obtained at long wave- 
lengths (A ^ £). For excitations where A is compara- 
ble to the size of the sample, this latter condition de- 
fines the Thomas-Fermi regime. Experiments have been 
performed either in the limit of zeroth sound excita- 
tions 1^,0, or intermediate to the two limits [§fi|. 

At higher thermal densities, when A 3> Z m f p collective 
excitations become hydrodynamic in nature, and one ex- 
pects two phonon-like excitations known as first and sec- 
ond sound. The presence of two hydrodynamic modes 
is similar to the case of superfluid 4 He. However, while 
in superfluid 4 He second sound is an entropy wave with 
no density oscillation, in a dilute Bose gas both first and 
second sound describe density waves. First sound, which 



is defined as having the larger speed of sound, involves 
mainly the thermal cloud, whereas second sound is sim- 
ilar in character to zeroth sound and is confined mainly 
to the condensate Jll],|l2j ■ 

In this Letter, we study collective excitations at non- 
zero temperature, and extend earlier work |(| in several 
ways. First, we study zeroth sound condensate oscilla- 
tions entirely in the Thomas-Fermi regime. In contrast, 
results of the earlier study were complicated by the tran- 
sition of the condensate from the Thomas-Fermi to the 
free-particle regime with increasing temperature. Sec- 
ond, near the critical temperature we approached the 
hydrodynamic limit and observed the onset of hydro- 
dynamic excitations analogous to first sound. Finally, 
we observed a new out-of-phase dipolar oscillation of the 
condensate and the thermal cloud, analogous to the out- 
of-phase second sound mode in liquid helium | jl2fl . 

The excitations were probed generally in three steps. 
First, as described in previous work p3| , we produced 
a magnetically confined, ultracold gas of atomic sodium 
in a well-determined equilibrium state. This state was 
controlled by adjusting the final frequency z/ r f used in 
the rf evaporation by which the clouds were cooled. The 
clouds were cigar-shaped with weak confinement along 
one direction (axial), and tight confinement in the other 
two (radial). Second, the cloud was manipulated with ei- 
ther time-dependent magnetic fields or off-resonant light 
to excite low-lying collective modes. Finally, the cloud 
was allowed to oscillate freely and probed in situ with 
repeated, non-destructive phase-contrast imaging [Q. 

To accurately characterize the magnetic trapping po- 
tential, we excited center-of-mass oscillations of the cloud 
in the axial direction by sinusoidally varying the position 
of the trap center. By observing the subsequent motion 
of the cloud, we measured the axial trapping frequency 
as v z = 16.93(2) Hz, while the radial frequency was esti- 
mated to be 230 Hz based on earlier measurements. 

The m = quadrupolar modes of the condensate and 
the thermal cloud were excited by a 5-cycle pulsed 
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modulation of the axial magnetic field curvature. Con- 
densate oscillations were produced with driving frequen- 
cies of 25 - 27 Hz, while the thermal cloud was excited 
at about 30 Hz. As shown in Figure 1, the most promi- 
nent feature of the condensate oscillations was the change 
in the axial length. Data were evaluated for oscillations 
with a relative amplitude of about 10% | Q |. 

The oscillations of each cloud were probed with 22 non- 
destructive images: one before the excitation to charac- 
terize the initial conditions, and three groups of 7 images 
during the oscillation. The three groups were separated 
by a delay time which was varied between 1 and 200 ms 
according to the damping time of the oscillation. This 
method of probing gave highly accurate single-shot mea- 
surements of oscillation frequencies and damping rates, 
allowed data to be collected efficiently, and overcame the 
additional fluctuations introduced when combining data 
from observations on several clouds. 

Phase-contrast images were analyzed by fitting the ob- 
served column densities h(r,z) 15 with the function 
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Here, r and z are the radial and axial coordinates, respec- 
tively, while l r and l z are the lengths of the condensate 
and ay and a z are the rms diameters of the thermal cloud. 
The function g^ is defined by gi[x) = X^i x V* 2 - Eq. 
|l| is motivated by a simple model of the mixed cloud: a 
mean-field dominated condensate amid a saturated non- 
interacting non-condensed Bose gas. Yet, since all quan- 
tities were allowed to vary independently, Eq. [I] is an 
almost model-independent parametrization of a bimodal 
distribution: l z is determined from the cusps of the bi- 
modal distribution, and ay from the thermal tails. 

The initial conditions for the oscillation were char- 
acterized by the total number of atoms N, the tem- 
perature T, and the chemical potential [i. N was ob- 
tained by integrating the column density, while T and fi 
were determined from the fits by fc^T = ir 2 mv 2 (j 2 and 
fi = ft mv 2 l 2 , where ks is Boltzmann's constant. For 
mixed clouds, T was determined by fitting the thermal 
wings alone. The initial conditions (Fig. 2) were varied 
by adjusting Aiyf = frf — Vhot where ^ ot is the resonant 
rf frequency at the bottom of the magnetic trap which re- 
mained constant within 20 kHz. The Bosc-Einstein con- 
densation transition was observed at T — 1.7 /iK with 
about 80 x 10 6 atoms. The temperature varied linearly 
with Af r f , with a slope of 3.5 /LtK/MHz, and was mea- 
surable down to 0.5 /iK. Reliable determinations of l z 
(and n) were obtained only below Aiy-f = 350 kHz. 
For low Ai/ r f , fi/ ks — 380 nK, corresponding to a con- 
densate of 15 x 10 6 atoms with a maximum density of 
no = 3.8 x 10 14 cm~ 3 using fi = Uq ■ 4nh 2 a/m with the 



scattering length a — 2.75 nm. 

The lengths l z and a z were fit independently to de- 
caying sinusoidal functions of time. We thus determined 
the frequency and damping rate of two distinct oscilla- 
tions. We observed no evidence for coupling between the 
oscillations in l z and in ay, allowing one to consider the 
excitations as nearly isolated oscillations of the thermal 
cloud and of the condensate, respectively. 

The thermal cloud oscillated at a frequency of about 
1.75 v z with a damping rate of about 20 s _1 , both above 
and below the transition temperature (Fig. 3). The ob- 
served frequency v is between the predictions of v = 2v z 
in the collisionless limit, and of v = 1.55 v z in the hy- 
drodynamic limit fl6|| . The damping rate is predicted to 
vanish in both the collisionless and hydrodynamic lim- 
its, and to reach a broad maximum of ~ 1.4 v when v is 
between its two limiting values [p"7f , in agreement with 
our observation. The collisional mean-free path Z m f p ~ 
/im x (T/fxK) 3 ' 2 using the peak density 
of the thermal cloud n T = 2.612 (mk B T/2Trh 2 ) 3 / 2 , and a 
collisional cross section a — 87ra 2 c . Around the transition 
temperature, we find <j z ~ 8 Z m f p . This comparison of 
length scales, the observed frequency shift away from 2iy , 
and the high damping rate all demonstrate that the col- 
lective behavior of the thermal cloud is strongly affected 
by collisions. Thus, the thermal oscillations which we 
observe indicate the onset of hydrodynamic excitations, 
i.e. first sound. However, the pure hydrodynamic limit, 
characterized by low damping, would only be reached for 
even larger clouds. 

We studied the quadrupolar condensate oscillations in 
greater detail. Typical oscillation data (Fig. 4.) demon- 
strate that the oscillation has a slightly lower frequency 
and is damped more rapidly at high temperature than 
at low temperature. At low-temperatures, the conden- 
sate oscillation frequency approached a limiting value of 
1.569(4) v z , close to the zero-temperature, high-density 
prediction of 1.580 v z ||. The slight difference between 
these values may be due to non-zero temperature effects 
even for small A^ r f . At higher temperatures, the fre- 
quency drops below the low-temperature limit. This 
trend might be explained by the simple non-zero tem- 
perature application of Bogoliubov theory [lq| , in which 
the condensate oscillates in a combination of the exter- 
nal trapping potential U(r, z) and a mean-field poten- 
tial exerted by the (static) thermal cloud. Considering 
a thermal density ut <x e - u ( r ^ z )/ k sT ^ one can es ti ma te 
that the mean-field potential shifts the collective excita- 
tion frequencies downward by as much as 5%. This shift, 
which scales as a/AdB where AdB is the thermal deBroglie 
wavelength, is consistent in magnitude with those which 
we observed. 

Damping rates for the condensate oscillations varied 
strongly with temperature, rising from a low-temperature 
limit of about 2 s^ 1 to as much as 15 - 20 s^ 1 at 
high temperatures. Recent treatments based on Landau 
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damping |19|-|21| provide qualitative agreement with our 
findings. A quantitative prediction for these damping 
rates |2l[] could not be checked because our data were 
collected for ksT < 6/z, where the high-temperature pre- 
diction might not be applicable, while the inability to 
measure T for low Ai/ r f prevents a comparison at low 
temperatures. 

Can the condensate oscillations be considered to be 
second sound? A comparison of length scales — l z — 
4/ m fp at high- temperatures — suggests that hydrody- 
namic effects may already be present. However, there are 
no theoretical predictions regarding the transition from 
zeroth to second sound with which to compare our data. 
In future experiments with larger condensates, the signa- 
ture of this cross-over may appear in the damping rate 
of the oscillations, which should decrease again at high- 
temperatures as one reaches the hydrodynamic limit. 

To further probe the interaction between the thermal 
cloud and the condensate, we studied an excitation of 
a different symmetry: the rigid-body, out-of-phase mo- 
tion of the condensate and the thermal cloud in the har- 
monic trapping potential JT^ ]. At temperatures where 
the condensate fraction approaches 50%, both compo- 
nents should participate equally in this motion. This 
mode is analogous to second sound in liquid helium, 
where the superfluid and the normal fluid undergo out- 
of-phase density oscillations of equal magnitude. 

We excited this mode by using off- resonant, blue- 
detuned laser light, which repels atoms due to the AC 
Stark shift fp|j7p^]. We produced a 3 /^K-high repulsive 
potential by focusing 40 mW of light at 514 nm to an 
elongated spot with 1/e 2 half- widths of about 12 /im and 
100 fim. After a partly condensed cloud was formed, the 
light was turned on and directed at the edge of the cloud, 
where it overlapped only with the thermal cloud. By tilt- 
ing a motorized mirror, the laser beam was steered to- 
ward and then away from the center of the cloud, thereby 
pushing the thermal cloud in the axial direction while 
not directly affecting the condensate. The light was then 
turned off, and the cloud allowed to freely oscillate. 

The position of the center-of-mass of the condensate 
(Fig. 5a) was monitored with repeated phase-contrast im- 
ages, and data from several repetitions of the experiment 
were collected with a variable delay between excitation 
and probing. The condensate motion was initially slow, 
and then grew to an asymptotic sinusoidal oscillation, 
corresponding to the in-phase motion of the entire cloud 
(condensate plus thermal cloud) in the magnetic trap at 
the trapping frequency v z . By subtracting the undamped 
center-of-mass motion of the entire cloud, we isolated the 
antisymmetric oscillation of the condensate and the ther- 
mal cloud (Fig. 5b). 

The frequency of the asymmetric dipole mode of 
17.26(9) Hz was significantly lower than the trapping fre- 
quency which was v z = 18.04(1) Hz at that time. This 
~ 5% frequency shift is again clear evidence of the in- 



teraction between the thermal cloud and the conden- 
sate. The motion exhibited in this oscillation can be 
regarded as the bulk flow of a superfluid through the 
surrounding normal fluid. Therefore, the damping which 
we observe directly measures the dissipation of superfluid 
flow. By varying the amplitude of the oscillation, one can 
study this dissipation as a function of velocity, and ob- 
tain more detailed information about the superfluidity of 
Bose-Einstein condensed gases. Further, the motion we 
observed, in which the condensate is driven by the mov- 
ing thermal cloud, cannot be described by recent theo- 
ries which assume a stationary thermal cloud ]l8[], and 
requires more sophisticated treatments p3] . 

In conclusion, we have studied the collective excita- 
tions of a dilute Bose gas at non-zero temperatures in the 
Thomas-Fermi limit, and near the hydrodynamic regime. 
The hydrodynamic oscillation of the thermal cloud, cor- 
responding to first sound, was indicated by measure- 
ments both above and below the Bose-Einstein conden- 
sation transition. The accurately determined frequency 
shift of the quadrupole oscillations away from their zero- 
temperature limit and of the antisymmetric oscillation 
away from the trap frequency are measures of the forces 
exerted by the condensate and the thermal cloud on one 
another. The damping of these modes sheds light on the 
dissipation of superfluid flow. 
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FIG. 1. In situ images of the m = quadrupolar con- 
densate oscillation. A Bose-Einstein condensate with no 
discernible thermal component was imaged every 5 ms by 
phase-contrast imaging. The evident change in the axial 
length of the condensate was used to characterize the os- 
cillation. Final data were evaluated for smaller oscillation 
amplitudes. 
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FIG. 2. Characterization of the trajectory across the 
Bose-Einstein condensation transition. The total number N 
(a) was determined by summing over the observed column 
densities. The temperature T (b, open circles) and chemical 
potential n (b, closed circles) were determined from fits to 
the data. These are plotted against Av r {. The dashed line 
indicates the observed transition temperature. 
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FIG. 3. Temperature dependent frequency and damping 
rates of m = quadrupolar collective modes. Points show 
measurements for oscillations of the thermal cloud (open 
circles) and condensate (closed circles). Horizontal dot- 
ted lines indicate the free-particle limit of 2 v z and the 
zero-temperature condensate oscillation limit of 1.580 v z . The 
vertical dashed line marks the observed transition tempera- 
ture. 
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a) Avrf = 30 kHz 
T<0.5|iK 
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Time After Drive (ms) 
FIG. 4. Damped quadrupolar condensate oscillations at 
low (a) and high (b) temperature. Points show the axial 
condensate length determined from fits to phase-contrast im- 
ages. Lines are fits to a damped sinusoidal oscillation with 
a downward slope. This slope accounted for heating during 
the oscillation. The oscillation at high temperature has a 
slightly lower frequency, and is damped more rapidly than at 
low temperature. 
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FIG. 5. Isolating the antisymmetric dipole mode. The cen- 
ter-of-mass of the condensate in the trap (a) was monitored 
after the selective displacement of the thermal cloud with a 
far-detuned laser beam. The oscillation of the condensate 
relative to the center-of-mass of the entire sample (b) oc- 
curred at a frequency below the trap frequency with a damp- 
ing rate of 9(l)s _1 . The initial conditions were T = 1/iK , 
fi/k B ^ 200 nK, and TV = 40 x 10 6 . 
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